NASA-CR-200728 


^7 


The Astrophysical Journal, 453; 464-479, 1995 November 1 

C 1995. The American Astronomicai Society. AJl rigbti reserved. Printed in U.S.A. 


SIMULTANEOUS EXTREME-ULTRAVIOLET EXPLORER AND OPTICAL OBSERVATIONS OF 

AD LEONIS: EVIDENCE FOR LARGE CORONAL LOOPS AND THE NEUPERT EFFECT IN 

STELLAR FLARES 

Suzanne L. Hawley* 

Dq^artment of Physics and Astronomy, Michigan State University, East Lansing, MI 48824 

George H. Fisher 

space Sciences Laboratory, University of California, Berkeley, CA 94720 

Theodore Simon 

Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

Scott L. Cully 

Space Sciences Laboratory, University of California, Berkeley, CA 94720 

SUSANA E. Deustua 

Institute of Geophysics and Planetary Physics, Lawrence Livermore National Laboratory, Livermore, CA 9455 1 

Marek Jablonski 

Institute of Theoretical Astrophysics, University of Oslo, N-031 5 Oslo, Norway 

Christopher M, Johns-Krull 

McDonald Observatory, University of Texas, Austin, TX 78712 

Bjorn R. Pettersen 

Geodetic Institute, Norwegian Mapping Authority, N-3500 Honefoss, Norway 

Verne Smith ^ 

Department of Astronomy, University of Texas. Austin, T7C 78712 

William J. Spiesman 

McDonald Observatory, University of Texas, Austin, TX 78712 
AND 

Jeffrey Valenti 

Joint Institute for Laboratory Astrophysics, University of Colorado, Boulder, CO 80309 
Received 1995 March 10: accepted 1995 May 12 

ABSTRACT 

We report on the first simultaneous Extreme-Ultraviolet Explorer {EUVE) and optical observations of flares 
on the dMe flare star AD Leonis. The data show the following features: 

1. Two flares (one large and one of moderate size) of several hours duration were observed in the EUV 
wavelength range; 

2. Flare emission observed in the optical precedes the emission seen with EUVE; 

3. Several diminutions (DIMs) in the optical continuum were observed during the period of optical flare 
activity. 

To interpret these data, we develop a technique for deriving the coronal loop length from the observed rise 
and decay behavior of the EUV flare. The technique is generally applicable to existing and future coronal 
observations of stellar flares. We also determine the pressure, column depth, emission measure, loop cross- 
sectional area, and peak thermal energy during the two EUV flares, and the temperature, area coverage, and 
energy of the optical continuum emission. 

When the optical and coronal data arc combined, we find convincing evidence of a stellar “Neupert effect” 
which is a strong signature of chromospheric evaporation models. We then argue that the known spatial 
correlation of white-light emission with hard X-ray emission in solar flares, and the identification of the hard 
X-ray emission with nonthermal bremsstrahlung produced by accelerated electrons, provides evidence that 
flare heating on dMe stars is produced by the same electron precipitation mechanism that is inferred to occur 
on the Sun. 

We provide a thorough picture of the physical processes that are operative during the largest EUV flare, 
compare and contrast this picture with the canonical solar flare model, and conclude that the coronal loop 
length may be the most important factor in determining the flare rise time and energetics. 

Subject headings: stars: chromospheres — stars: flare — stars: individual (AD Leonis) — ultraviolet: stars 
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1. INTRODUCTION 

AD Leonis is a dMe flare star located in the immediate solar 
neighborhood, at a distance of about 5 pc. According to Petter- 
sen, Coleman, & Evans (1984), detectable flares (Johnson 
C/-band enhancements of ^0.1 mag) occur every 100 minutes 
or so; more energetic flares occur with an exponential decrease 
in frequency. Since AD Leo is the brightest active flare star 
visible from the northern hemisphere, it has been the target of 
several previous multiwaveicogth observing campaigns (e.g., 
Rodono et al. 1989; Hawicy & Pettersen 1991 [hereafter 
HP91]; Bookbinder, Walter, & Brown 1992). We report here 
on a successful observing campaign using the NASA Extreme- 
Ultraviolet Explorer (£C/K£) satellite and six ground-based 
optical telescopes. To our knowledge these are the first simul- 
taneous optical and EUV observations of stellar flares with the 
£t/F£ satellite. 

The effort involved in orgamzing such an extensive obser- 
vational campaign is more than justified by the importance of 
observing flare emission from both the high-temperature 
(coronal. X-ray) and low-temperature (chromospheric, optical) 
portions of the flaring plasma simultaneously. Previous inves- 
tigations of this type arc few and have reported varying results. 
Some workers have found that the coronal and chromospheric 
plasmas respond on a similar timescale to flare heating (e.g., 
Kahler et al. 1982; Haisch ct aL 1983; de Jager et al. 1989), 
leading to the relations of Butler, Rodono, & Foing (1988), 
Haisch (1989), and HP91 between the energy emitted in high- 
temperature (X-rays) and low-temperature (hydrogen and 
calcium lines, f/-band continuum) flare features. However, 
other workers (Doyle et al. 1988; Haisch et al. 1981) found that 
there was no simultaneous response between the high- and 
low-temperature plasma during the flares that they observed. 

Our goal was to obtain a new set of simultaneous coronal 
and chromospheric flare observations with higher spectral and 
temporal resolution than was previously attainable, and to 
analyze the data within the context of the generally accepted 
model for solar flares (e.g., Lin & Hudson 1976; Dennis & 
Schwartz 1989). In this modd, a flare begins with the rapid 
acceleration of a burst of nontfaermal electrons down the mag- 
netic field lines of a closed coronal loop (or an arcade of loops). 
The electrons collide with ions in the ambient plasma and emit 
hard X-rays via nonthermal brcmsstrahlung, from which a 
power-law distribution can be inferred for the electron energy 
spectrum. White-light continuum emission is often observed 
from small kernels at the footpoints of flare loops. This emis- 
sion is temporally and spatially correlated with the hard X-ray 
emission (see, e.g.. Fig. 3.25 of Canfield et al. 1986; Hudson et 
al. 1992) and is therefore emitted at or near the sites where the 
accelerated electrons impact the preflare chromosphere and/or 
photosphere (although the exact emission mechanism remains 
unknown). T^e spiraling motion of the energetic electrons 
around the magnetic field lines also results in gyrosynchrotron 
emission at microwave frequencies. The nonthermal electrons 
heat the lower atmosphere, resulting in both emission from the 
low-temperature plasma and the “evaporation’' of chromo- 
spheric plasma to coronal temperatures, filling the coronal 
loop with hot plasma. The term “evaporation ” has been tradi- 
tionally used to describe the “phase change” that occurs when 
dense, cool chromospheric plasma is heated during the flare 
and is transformed into hot, tenuous coronal plasma. As the 
corona increases in density (column depth), it b^omes increas- 
ingly effective at stopping the energetic electrons, which now 
heat the corona directly, while the amount of heat deposited in 


the lower atmosphere diminishes. During this time the corona 
emits thermal soft X-rays. When the coronal heating ceases, 
the amount of material in the loop reaches its peak value and 
then gradually condenses back to its preflare configuration, 
while the soft X-ray flare emission dies away. 

This canonical solar flare model manifests itself empirically 
through the well-known “Neupert effect” (Neupert 1968; 
Dennis & Zarro 1993). Because chromospheric evaporation, 
which is responsible for increasing the coronal emission 
measure, is proportional to the flare heating rate, which is in 
turn proportional to the hard X-ray flux, one anticipates that a 
time integral of the hard X-ray flux should be proportional to 
the observed thermal soft X-ray flux. Dennis & Zarro (1993) 
examined a data set of flares seen in both hard X-rays with the 
HXRBS spectrometer on SMM and in soft X-rays with the 
GOES satellite. They found that more than 80% of large events 
exhibit the Neupert effect. In this paper, we will use our optical 
data (as a proxy for hard X-rays)^ and EUVE (essentially soft 
X-ray) data to search for a stellar flare analog to the Neupert 
effect. 

In previous work we have quantified the overall theoretical 
framework of the canonical solar flare model by explicitly 
solving the energy equation for loops in the evaporating and 
condensing regimes (Fisher & Hawley 1990, hereafter FH90). 
Given a prescribed flare heating rate, the loop evolution can 
then be determined. In this paper, we extend the analysis in 
FH90 so that the new EUVE data can be used t9 cietermine 
parameters such as the loop length, maximum column depth (a 
measure of the amount of evaporated material), and coronal 
pressure, as well as the flare area coverage, emission measure, 
and toti thermal energy. We also use the techniques devel- 
oped in Hawley & Fisher (1992, hereafter HF92) to analyze the 
optical flare emission and to determine the temperature, area 
coverage, and total energy deposited in the lower atmosphere. 
The optical and coronal results are then combined to investi- 
gate properties of the flare heating mechanism and the efficacy 
of the solar flare picture in describing the observed stellar 
flares. 

2. THE DATA 

As part of the first year of the EUVE guest investigator 
program, we obtained EUV observations of AD Leo over 42 
satellite orbits between 1993 March 1-3. The data consist 
of spectra in the short-, medium-, and long-wavelength regions, 
plus EUV “photometry ” obtained with the Deep Survey (DS) 
instrument. We defer discussion of the spectra to a later paper 
(Cully et al 1995) and concentrate here on the DS data. To 
complement the EUV observations and provide information 
on the chromospheric flare response, we organized a simulta- 
neous optical campaign, which resulted in data from six optical 
telescopes: the Lick Observatory 24 inch (0.61 m), 

(photometry), 40 inch (1.02 m) (spectroscopy), and 24 inch (0.61 
m) Coude Auxiliary Telescope (echelle spectroscopy); the 
McDonald Observatory 36 inch (0.92 m) (photometry) and 82 
inch (2.08 m) (echelle spectroscopy); and the 100 inch (2.54 m) 
Nordic Optical Telescope (spectroscopy). 

^ With existing telescopes, hard X-ray emission is unobservable from stellar 
flares, even for the nearest stars, because nonthermal brcmsstrahlung is such an 
inefficient process. While impulsive gyrosynchrotron emission from highly 
energetic electrons may have been detected in a few stellar flares (see § 4 3 of 
the review by Bastian 1990; Fig. 1 of Rodono et al. 1989), the white-light 
continuum is the only easily measurable tracer of the impulsive energy release 
on flare stars, if we adopt the solar model. 
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2.1. Observational Ensemble 

Table 1 is an observing log for the seven telescope/ 
instrument combinations including the observers, type of data 
obtained, and observing intervals. We discuss each individ- 
ually below. 

2 . 1 . 1 . EVVE/DS 

Malina & Bowyer (1991) provide a description of the EUVE. 
Briefly, the telescope has a 2?1 diameter field of view and an 
on-axis angular resolution of ^ 18^. A specific target is visible 
for some 30 minutes of each 90 minute orbit, with visibility 
defined by Earth occultation and the geocoronal background 
count rate. The Deep Survey (DS) instrument is also described 
in Malina & Bowyer (1991) and in Siegmund et al. (1986). We 
used the Lexan/Boron filter, whidi has a bandpass of 65-190 A 
at 10% of peak transmission. The spectroscopic data (Cully et 
al. 1995) show that the strongest lines in this bandpass are 
mainly from highly ionized stages of Fe formed at tem- 
peratures of 10^ to several x 10^ K. 

The DS data were processed into files containing all counts 
detected within 12' of the target portion, along with the arrival 
time, detector location, and aspect-corrected sky position of 
each photon. The average source count rate was higher than 
the local sky and detector background rates, so the source 
location could be determined visually. We determined that the 
spacecraft pointing during our AD Leo observation was stable 
to within 5" over more than 90% of our exposure time, and the 
occasional short deviations (all <70*) that did occur took 
place when AD Leo was in a state of low activity. A sym- 
metrized version of the DS point spread function was used to 
calculate the source and background count rates. The data 
were then corrected for instrumental deadtime and 
“ Primbsch ” effects (Abbott, Miller, & Gavin 1993) and binned 
into 500 s intervals. The optically bright binary giant, y Leo 
(KO III -I- G7 III), lies 6' to the SE of AD Leo. We have exam- 
ined a deep, 26 ks exposure of AD Leo taken with the PSPC 
instrument aboard ROSAT (image rp200076) which we 
obtained from the ROSAT archives, and we find that neither of 
the late-type giant stars is present as a soft X-ray source. 
Therefore we attribute all the EUVE emission we have 
observed to AD Leo. 

A possible source of uncertainty in the DS count rate is 
introduced by an area of decreased sensitivity in the center of 
the DS microchannel plate detector that appeared after an 


extensive observation of the white dwarf HZ 43 in 1993 Feb- 
ruary (Miller 1994). This **dead spot” is centered on the bore 
site of the instrument, with a maximum loss of efficiency in 
some channels of ^80%. To determine whether the dead spot 
affected the DS count rate in our observation, we compared 
the wavelength-integrated count rate in the SW detector of the 
EUVE spectrograph with the DS count rate and found that 
the count rates are consistent with one another if no dead spot 
coirection is applied to the DS data. We also compared the 
ratio of the DS count rate to the SW count rate from the 
observation <rf AU Mic (made before the observation of HZ 43) 
(Cully et aL 1993, 1994) to the same ratio from our AD Leo 
obsemtions and found that it did not change, further indicat- 
ing that the DS count rate was not strongly affected by the 
dead spot during our AD Leo observations. Therefore we have 
not applied any dead spot correction to our DS count rates for 
AD Leo. If it becomes clear at a later date that the count rates 
we report here need to be changed by a “ correction factor 
then some of the coronal parameters that we derive in § 4.2 for 
the flare observations will be affected. The footnotes to Tables 
2 and 3 describe how the quoted values scale with The 
valued, = 1 adopted here is probably known only to within a 
factor of 2 or so based on the tests we have performed. Values 
of Alt as large as 3-4 (i.e., corresponding to the possible 
maximum loss of efficiency of ^ 80%), while changing some of 
the coronal parameters, will not significantly affect our conclu- 
sions. ' ' 

2 . 1 . 2 . Photometry 

Photometric observations were made at the Lick Obser- 
vatory 24 inch telescope and the McDonald Observatory 36 
inch telescope. The photometry at Lick was obtained with a 
single-channel photometer, which was equipped with standard 
Johnson C/, B, and V filters and a 1P21 photomultiplier cooled 
with dry ice. Because the filters had to be changed manually, 
AD Leo was monitored only in f/, with occasional obser- 
vations in B and V. The integration time was 2 s in each filter. 
The observing periods were photometric, and the data were 
corrected for extinction and reduced to the standard system 
using observations of standard stars (Landolt 1973, 1992). 

At McDonald, a two-star photometer was used for UBV 
photometry. The primary photomultiplier was thermoelcctri- 
cally cooled while the secondary tube was operated at ambient 
temperature. AD Leo was observed through the primary 


TABLE 1 


Observing Loo 


Telescope/ 

Observm 

Instrument 

Wavdengtbs/ 

Resolution 

Tune (UT) 

Elapsed 
Time (hr) 

EVVE/Simon 

Deep Survey 

65-190 A 

Mar 1 6.0- 
Mar 3 23.3 

6.0-71.3 

Lick 24 inch/Johos, Hawley 

photometer 

U{BV) 

Mar 1 7.S-9.5 
Mar 2 5.0-12.1 

7.8-9.5 

29.0-36.1 

Lick 40 inch/Deustua, Hawley 

medium-resolution 

spectrograph 

3900-7000 A 
S A 

Mar 2 7.4-11.7 

31.4-35.7 

Uck CAT/CuUy, Valenti 

Hamilton echelle 

4800-0700A 

Mar 1 5,0-110 

5.0-110 


spectrograph 

R - 4S.000 

Mar 2 3.5-115 

27.5-36.5 

McDonald 36 inefa/Spiesman 

photometer 

UBV 

Mar 2 8.3-10.3 
Mar 3 5.9-10.5 

313-34.3 

53.9-58.5 

McDonald 82 inch/Smith 

Sandiford echelle 

600O-M0OA 

Mar 2 9.7, 11.3 

33.7, 35.3 

NOT 100 inch/ 

spectrograph 

K = 60.000 

Mar 3 6.2-7.0 

54.2-55.0 

Pettersen, Jablonski 

medium-resolution 

3900-4900 A 

Mar 2 0.7-4.6 

24.7-28.6 


spectrograph 

3.8 A 

Mar 2/3 21.5-1.4 
Mar 3 S.9-6.6 

45.5-49.4 

53.9-54.6 
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channel, and the second channel was used to monitor sky con- 
ditions. An integration time of 2 s was used at each of the f/, £, 
and K filter positions. An automatic filter wheel and computer- 
controlled integrations provided a total cycle time of 7 s. The 
sky was photometric for most of the observing period, 
although high winds forced occasional temporary closures of 
the dome. As with the Lick data, Landolt standards were used 
to correct for atmospheric extinction and to reduce the data to 
the standard UBV system. The absolute photometric measure- 
ments at the two sites agreed to within -^0.02 mag in each filter 
during periods of simultaneous observing, indicating external 
agreement at the level of 2%, 

2. 1 .3. Medium-Resolution Spectroscopy 

Medium-resolution spectroscopy was obtained with the 
Lick 40 inch Nickel Telescope and with the 100 inch Nordic 
Optical Telescope (NOT) at La Palma. The Lick observations 
were made with a Cassegrain Spectrograph equipped with a 
600 lines mm * ‘ grating, giving an effective resolution of 2.6 A 
pixeP ‘ on a 400 x 1200 Reticon CCD detector. The slit width 
was 3", corresponding to 2 pixels on the CCD. The spectra 
cover the wavelength region 3900-7000 A. Exposure times 
varied from 20 to 30 s, with a minute between successive expo- 
sures for readout and recording Calibration frames included 
an argon-neon arc lamp for wavelength calibration and spectra 
of standard stars from Massey et aL (1988) for flux calibration. 
The data were reduced and transformed to a calibrated flux 
system using IRAF. We expect that there were small (;S5%) 
light losses from the edges of the slit as a result of seeing 
changes and guiding errors. The night was photometric. 

The observations on the NOT were made with the Low 
Dispersion Spectrograph. The spectra cover the wavelength 
region from 3900 to 4900 A with an effective resolution of 1.9 A 
pixel" ^ The slit width was ITl which corresponded to 2 pixels 
on a thinned, back-illuminated Tektronix CCD. The exposure 
time was 20 s, with 5 s between successive exposures for 
readout and recording. The data were reduced using routines 
written in the Interactive Data Language (IDL). The weather 
was not photometric, so only the equivalent widths of the prin- 
cipal emission lines were obtained. 

2.1.4. Echelle Spectroscopy 

High-resolution spectra were obtained in the 4800-8700 A 
region at Lick Observatory with the 24 inch Coude Auxiliary 
Telescope (CAT) and in the 6000-8000 A region at McDonald 
Observatory with the 82 inch Struve Telescope. The CAT 
observations employed the Hamilton Echelle Spectrometer 
and a cooled TI 800 x 800 CCD (Vogt 1987). The spectro- 
meter was designed for use with a 2048 x 2048 CCD; thus, use 
of the TI 800 x 800 CD resulted in gaps in the spectral cover- 
age. The resolving power of the spectrometer is approximately 
R = X/A2 = 48,000; there were 2.5 pixels per resolution 
element. The slit width of 500 ftm that was used corresponds to 
1.5 pixels on the detector and 477 on the sky. The data were 
reduced using the IDL procedures outlined by Basri, Wilcots, 
& Stout (1989) and Valenti (1994). Exposure times were 15 
minutes per frame during quiescence and 5 minutes per frame 
during flare. The readout time was approximately 1 minute for 
each exposure. 

Observations at McDonald were made with the Sandiford 
cross-dispersed echelle spectrometer (McCarthy et ai. 1993) at 
the Cassegrain focus. The detector is a Reticon 400 x 1200 
CCD with 27 pixels; the two-pixel resolution is 60,000. The 
two-dimensional raw frames were reduced to single-order one- 


dimensional spectra using the IRAF echelle data reduction 
routines. Integration times were 2 minutes, with each extracted 
spectnun having a signal-to-noisc ratio of about 40 in the 
** continuum ” at Ha. 

2.2. EUV and Optical Flare Observations 

Figure 1 shows the EUV light curve obtained with the 
£(/K£/DS instrument. Each point represents an independent 
bin of 500 s of integration time. The time axis is measured in 
hours from 1993 March 1.0 UT. Throughout the paper, we will 
use this time axis as a reference for all the data, including the 
optical (see col. [5] of Table 1). Also indicated in the lower half 
(^Figure 1 are the observing intervals for the optical data. It is 
dear that two major EUV flares occurred during the time we 
were monitoring the star. The first (hereinafter EUV flare 1, or 
EFl) began between 33:30 and 34:30 hr, reached a peak count 
rate of ^1.0 counts s" \ and continued for at least 8 hr, while 
the second (hereinafter EUV flare 2, or EF2) began between 
55:30 and 56:30 hr, reached a peak count rate of ^0.8 counts 
and continued for approximately 3 hr. Note that EF2 
occurred while the DS count rate was still elevated from the 
decay of EFl. We were fortunate to have optical coverage 
which overlapped both flares, during the time intervals 
denoted in Figure 1 by C (EFl) and E (EF2). TTiere are three 
other intervals with simultaneous optical and EUV obser- 
vations, and these are labeled A, B, and D in the same figure. 

Figures 2a-2e show the data obtained during the five simul- 
taneous observing periods, A-E. Period A is a time of relative 
quiescence, while B and D show small flares with little or 
no EUV response. Note that we are unable to resolve low-level 
flare enhancements in the EUV data because of the possible 
small-scale variability introduced by the presence of the dead 
spot. Thus, we can make no quantitative statements about the 
presence or absence of low-level, short-lived EUV emission 
such as might be expected to accompany these small optical 
flares. 

The first large EUVE flare, EFl, begins during period C. 
Both photometric sites observed a (/-band flare beginning at 
32:58 hr, which reached a peak brightness some 0.7 mag above 
the quiescent level, followed by a sharp and then more gradual 
decline. The (/-band flux did not return to its quiescent level, 
but remained enhanced for several hours. A second (/-band 
flare occurred at 35:48 hr, peaking 0.6 mag above the preflare 
level, or equivalently 0.8 mag above the quiescent level. Our 
medium- resolution spectroscopy shows that the hydrogen 
lines followed the (/-band flare in the expected way, increasing 
in flux during the onset of the initial flare, followed by a slow 
decline. Figure 2c (bottom panel) gives the Hy light curve at a 
time resolution of 5 minutes. We also have data for the other 
Balmer lines, which follow the same time evolution as Hy, and 
for Ca II K, which is quite noisy but appears to have a some- 
what slower evolution, in accord with the results in HP91. The 
Balmer decrement at the peak of the first (/-band flare 
(preceding EFl) is consistent with the decrement obtained by 
HP91 and hence, as shown there, with other stellar flares. The 
flaie decrement is quite shallow compared to the quiescent 
decrement, indicating that the luminosity in the higher Balmer 
lines increases markedly compared to Ha. 

The Balmer decrement provides an observational constraint 
on chromospheric models (such as those in Houdebine & 
Panagi 1990; HF92; Houdebine & Doyle 1994a, b), particu- 
larly when information on the flare heating, through optical 
continuum and high-energy observations, is available. In solar 
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Fig. I . — EC/KE/DS light curve for the entire period of observation, with lime (in hr) measured from 1993 March 1.0 UT. Horizontal bars at the bottom of the plot 
indicate the times when coordinated observations were being nuuJe from each of the ground-based telescopes in the campaign. Periods of overlapping observations 
are labeled A-E in chronological order. 

HP equivalent width. The Ha data begin shortly after the flare 
peak but arc enhanced. Another If -band flaring episode begins 
at 55:25 hr and continues with short timescale variability for at 
least 3 hr until the end of the observations. Overlapping this 
optical activity, EF2 begins sometime between 55:30 and 
56:30 hr, reaches a peak that lasts at least 1.5 hr l^tween 56:30 
and 58:00, and decays by 60:00 hr. There is not a clearly 
measurable time lag between the optical and EUV emission in 
EF2, possibly because of our inability to resolve the EUV data 
on the shorter timescales that arc involved. 

In Table 2 we compile the empirical results from the obser- 
vations, including the rise and decay phase timescales, count 
rates, average magnitudes, and total (time integrated) energies 
for both the optical and EUV passbands in each flare. The 
EUV energy was computed from the count rate by assuming 
an average photon wavelength of 1(X) A and an effective 
EUVE/DS detector area of 25 cm^ (Cully et al. 1994, hereafter 
CFAS, Fig. 2), while the L/-band energy comes from the mea- 
sured magnitudes and the usual Johnson (1966) calibration. 
Note that since our C/-band data fail to cover the entire period 
of the EFl flare, our energy estimate is a lower limit to the total 
Lf-band energy radiated over the course of the flare event. Both 
the optical and EUV energy estimates assume that the flare 
radiation is emitted isotropically from a source that is 
unocculted. 

3. DISCUSSION OF THE OPTICAL DATA 

The optical continuum of AD Leo, as measured by the 
l/-band photometry, shows long-lived flaring episodes during 
both EUV flares. In most previous work on stellar flares, the 


flares, Ha becomes greatly enhanced. Stellar flares may thus 
differ from their solar counterparts in this regard. To determine 
whether such behavior reflects merely a difference in atmo- 
spheric structure, or a more fundamental change in the type 
and/or amount of chromospheric flare heating, will require 
more detailed modeling of the flare chromosphere. 

Unfortunately, the star set at Lick shortly after the second 
optical flare, and we have no further optical data for the dura- 
tion of EFl. 

The top panel in Figure 2c shows the EUVEfDS data during 
the overlapping time period. It is immediately apparent that 
EFl is not coincident with the beginning of the optical activity. 
In fact, the EUVE orbit spanning 32:50-33:20 hr exactly 
covers the onset of the first optical flare, and no simultaneous 
EUV flare was observed. The minimum time delay between the 
optical and EUV flares is ^25 minutes, corresponding to the 
end of this orbit, when EUV data were no longer being 
acquired. Clearly by the start of the next orbit (at 34:30 hr), the 
EUV count rate was enhanced, representing the first data 
obtained during the rise phase of EFl. The first noteworthy 
result of this work is, ther^ore, the considerable delay between 
the rise in optical light (both continuum and emission lines) 
and in the EUV photometry. A peak in the DS data is seen 
during the orbit near 36:15, and for the next several hours the 
flare emission hovers near this peak value. A sharp decay from 
the flare peak is clearly visible in the orbit near 42:00, and a 
change to a shallower decline occurs near 44:00, where the DS 
counts are still elevated above the quiescent value. 

Flare EF2 occurred during period E (Fig. 2e). Just after 
54:00 hr, an optical flare is seen in both the U-band and the 
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“canonical” l/-band light curve has a rapid rise of some 
seconds to a few minutes, which is followed by a slower, expo- 
nential decline lasting some tens of minutes (Moffett & Bopp 
1976). In general there is also a correlation between the rise 
and decay times and the flare energy, in the sense that more 
energetic flares take longer to rise and decay. The flares report- 
ed here do not follow such a neat picture. The two distinct 
optical flares observed during period C (accompanying EFl) 
do show rapid rise times, but the decay from the first peak is 
interrupted, remaining elevated at ^0.2 mag above the quies- 
cent level for almost 3 hr before rising to the second peak. The 
majority of the flare energy is emitted not during the peak, but 


during the long period of decline. Such behavior, though 
unusual, is by no means unprecedented. Long-lived, low-level 
flaring episodes have been reported previously by a number of 
authors (e.g., Byrne 1979, Fig. 3; Moffett 1974, Tables 17-23; 
Pcttersen et al. 1986, Table III and Fig. 5). Also, there has been 
a tendency in the past to ignore continuing low-level emission 
after the main flare peak has passed (see, e.g., Moffett & Bopp 
1976, Figs. 1-3), so that extended periods of activity may 
actually be more common than is generally thought. The 
advantage of having two photometric telescopes observing 
simultaneously during period C makes our observations 
unequivocal; the excess emission cannot be ascribed simply to 
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TABLE 2 


Measured Quantities 


Quantity 

EFl 

EF2 

AU Mic 


EUV 



Rise time (s) 

Decay time (s) 

Preflare count rate* 

c^(»-‘) 

Peak flare count rate* 

(s-‘) 

Flare decay count rate* 

C.(s-‘) 

Total energy 

Radiated £ds (ergs)* ** 

11000 

6000 

021 

0.98 

0.65 

4.1 X 10^* 

<6000 

<6000 

0.46 

0.76 

0.46; 

8.2 X 10^* 

5600 

4200 

0.50 

7.00 

4.30 

3.9 X 10“ 


Optical 



U magnitude, preflare 

U magnitude. 

IZOO 

1100 


flare average 

11.70 

11.82 


B magnitude, preflare 

B magnitude. 

10.98 

10.98 


flare average 

10.90 

10.93 


Flare time interval (hr) 

Total U flare 

319-36.1 

55.4-58.4 


energy, £*; (ergs) 

>4.9 X 10^* 

16 X 10“ 



• The count rates and DS energy listed here assume no correction for 
possible “dead spot ” effects in the AD Leo observation. If a correction factor 
/d, is applied to the DS count rate, each quantity should be multiplied by/^. 

** £ds is the total energy radiated during the rapid rise, flare peak, and rapid 
decay phases which we model in this paper. It does not include the ener^ 
radiated during the long slow decay phase seen after the rapid decay phase in 
EFl and the AU Mic flare. 

instrumental effects. In fact, careful optical monitoring might 
reveal that activity of this type is actually prevalent on dMe 
stars. The two optical flare events accompanying EFl and EF2 
each have total U-band energies of >10“ ergs (see Table 2); 
thus, they are more energetic than ~ 95% of the flares studied 
by Pettersen et al. (1984). 

The detailed 17-band flare light-curve morphology 
undoubtedly contains interesting information about the flare 
heating in the deep atmosphere. Unfortunately, until we can 
identify and model the exact mechanism by which the U-band 
continuum is produced, we are unable to use the detailed light- 
curve morphology to study flare heating in a quantitative way. 
We do, however, use the overall shape and energy to investi- 
gate the Neupert effect during stellar flares (see § 5). 

3.1. Diminutions in the Stellar Continuum (DI Ms) 

Figure 2c shows several other intriguing features in the con- 
tinuum behavior that arc corroborated by observations at 
more than one telescope. Diminutions in the stellar continuum 
appear at times 33 : 52, 34 : 1 5, and 35 : 27. These diminutions (or 
DIMS, as we choose to call them) may be similar to the “ dips ” 
first described by Cristaldi, Gershberg, & Rodono (1980), and 
also reported by Giampapa et al. (1982), Doyle ct al. (1988), 
and Rodono et al. (1989). Haisch et al. (1983) discussed several 
empirical explanations for these features. Grinin (1983) sug- 
gested that such features may be produced when the initial 
burst of flare energy impacts the photosphere, leading to a 
sudden rise in the H' opacity and momentarily upsetting the 
radiative equilibrium. Henoux et al. (1990) discussed possible 
“black light flares" on the Sun using a similar model. Gur- 
zadyan (1988) proposed a way to distinguish DIMs obscr- 
vationally by observing both the blue and red response— the 


near-IR flux should become brighter during the time that the 
Uue flux is diminish ed- Interestingly, the reverse of this pheno- 
menon (IR dips during U-band enhancements) is seen in some 
flare observations (Rodono et al. 1989). Geometric effects (e.g., 
prominences blocking part of the stellar surface) or even instru- 
mental proWems cannot be entirely ruled out. We are con- 
fident that the long-lasting diminutions at ~33:52 and 
~ 35:27 are real, but the short, deep excursion at ~ 34: 15 seen 
at both photometric sites might be coincidental, since the 
shape and timing of the event differ when the two sets of data 
are examined in detail. The fact that both longer lived diminu- 
tions occurred just prior to optical flares may also be impor- 
tant and fits with both previous observations and Grinin’s 
theory. 

3.2. Temperature, Area Coverage, and Energy 

In HF92 we analyzed the continuum emission from the very 
large 1985 April 12 flare observed on AD Leo and showed that 
the optical continuum observed in four bandpasses (U, B, V. 
and R), as well as the UV continuum seen in the SWP and 
LWP bandpasses on the International Ultraviolet Explorer 
(/UE), could be fitted by a blackbody spectrum at a tem- 
perature of ~9000 K, with only a small area coverage factor of 
the star. Other proposed continuum emission mechanisms. 
ifiHnHIng free-free emission and hydrogen recombination radi- 
ation, could not explain the observed continuuni spectrum. 
Several other authors (Mochnacki & Zirin 1980; Giampapa 
1983; Kahler et al. 1982) have also fitted a blackbody spectrum 
to optical flare continuum emission. The current data for AD 
Leo U and B observations during the optical flares (the 

V band is not significantly enhanced relative to the quiescent 
flux for these flares), which are insufficient to make an indepen- 
dent determination of the emission mechanism. Thus, we 
assume that the emission comes from a blackbody source, as in 
the previous flares, and we solve for the temperature and area 
coverage of the source. 

The quiescent energy fluxes observed at earth m the U and B 
bandpasses prior to a flare are designated by Ei,„ and Fg„. We 
ag<iitiM! that during the flare a fraction X of the stellar area 
emits as a blackbody at temperature TJ,. The flux radiated in 
the U and B bandpasses by the flaring areas of the star is then 

convolved with the respective filter response functions. 
If Et^TJi) and EafTf,) are the fluxes in the U and B bandpasses 
that would be observed at Earth if the entire star were at 
temperature T„ (the convolutions described above diluted by 
the factor Rj/d^), then the increases in the U and B fluxes AF^; 
and AFg observed at Earth during the flare can be written 

AFu = XFtATn)- XF^,, ( 3 . 1 ) 

and 

AEg = XEglTr,) - XFg, . (3.2) 

These equations do not include possible projection effects or 
partial occultation of the flaring area behind the limb, and so 
the derived area coverage, X, should be used with appropriate 
caution. 

A transcendental equation for the flare temperature /fi can 
be found by dividing the two equations, resulting in 

AFc; ' F^;(rn) - F^, ‘ 

Once 7J| has been determined from equation (3.3), the coverage 
factor X can be found from either equation (3.1) or (3.2). 
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The quiescent magnitudes for AD Leo (see Table 2) are 
f/o = 12.00 and Bq = 10.98. During the rise phase of EFl, the 
average flare magnitudes are 11.70 and Bn ^ 10.90. 
Equation (3.3) then gives 7Jj ^ 9000 K (compatible with esti- 
mates made for other flares) and X 1,0 x 10'’^ An estimate 
of the total flare optical luminosity is then given by L„p, = (tT^ 
X X X 4nRl ^ 4 X 10^’ ergs s"^ The flare emission con- 
tinues for at least 3.2 hr (when the observations ended), so a 
lower limit to the optical flare energy is ^ 4.6 x 10^^ ergs, 
or about 10 times the measured l/-band energy (Table 2). Flare 
EF2 has average flare magnitudes of l/n ^ 11.82 and Bn ^ 
10.93, giving 7„ - 9500 K and X - 5.3 x 10■^ and - 2.5 
x 10^® ergs s" *. In this case, we estimate an optical flare dura- 
tion of 3.1 hr, from which we find 2.8 x 10^^ ergs, again 
approximately 10 times the C/-band energy in Table 2. These 
derived optical flare quantities are given in Table 3, along with 
the coronal properties that we obtain from the models in the 
next section. 

4. EUV DATA AND CORONAL LOOP MODELS 

Using the EUV flare data presented in § 2, we now seek to 
apply elements of our flare loop evolution model to determine 
the physical parameters of the coronal loops within which 
flares EFl and EF2 occurred. First we review the basic pre- 
cepts of the model and derive expressions relating observed 
quantities (such as the flare rise and decay times) to the model 
parameters (such as the loop length). Then we analyze flares 
EFI and EF2 and discuss the resulting loop parameters with 
reference to previous work. 

4.1. Flare Loop Evolution 

In FH90, we solved an energy equation for a coronal loop in 
the evaporating and condensing regimes. In order to use an 
analytic expression for the energy balance in the loop, it was 
necessary to make several basic assumptions about the loop 
conditions and geometry. Our model is strictly valid only when 
these assumptions are realized. The assumptions are as 
follows: (1) The gas pressure is small compared to the magnetic 
pressure within the loop; (2) Flare timescales are long com- 
pared with the acoustic transit time along a coronal loop, so 

TABLE 3 

Derived QuANTmei 


Quamiiy EFi EF2 AU Mic 


Optical 


T„ (K) 9000 9500 

1.1 X 10** 5.6x10*'' 

t o X 10'* 53 X 10“ * 

(ergs) >4.6 X 10** IS x 10** 


Coronal 


L(cm) 3.8 x 10*® <1.5 x 10*® 2.6 x 10*® 

1.3 x 10** <Ij0x 10** 1.5 x 10** 

(dyne cm “ *) 180 <280 350 

9.1 X 10*® 1.9 x 10*® 1.7 x 10** 

0.0085 0n0l8 0.045 

V (cm*)* 7.1 x 10*® 5.6 X 10*® 8.8 x 10** 

(cm *)* 8.2 X 10** IS x 10*' 2.9 x 10** 

£,h (ergs)* 1.9 X 10** 23 x 10** 4.6 x 10*® 


* The coronal area coverage, volume, emission measure, and thermal 
energy listed here assume no correction for possible **dead spot effects in the 
AD Leo observation. If a correction factory^ is applied, each quantity should 
be multiplied by 


that the corona is regarded as always being in a state of quasi- 
steady hydrostatic equilibrium; (3) The pressure scale height 
Hp is lai^ compared to the loop semi-length L, so that the 
pressure in the loop is uniform, although changing with time; 
(4) The loop has a uniform cross-sectional area, which does not 
vary with time; (5) The conductive heat flux is given by the 
Spitzer formulation (Spitzer 1962); (6) Spatial variations of the 
flare heating rate in the corona can be averaged, i.e., the volu- 
metric flare heating rate is assumed to be uniform in space 
(although it can vary in time); and (7) Radiative cooling takes 
place in the optically thin “ coronal ” approximation, with the 
temperature dependence of the cooling rate described by a 
single power-law fit over the temperature range of interest (see, 
e.g., the fit oa p. 662 of Rosner, Tucker, & Vaiana 1978, here- 
after RTV, to the cooling law of Raymond & Smith 1977). 

Assumption (1) is true in the solar case and appears reason- 
able for stara l^t facto, we find that assumptions (2) and (3) 
are satisfied by the flare evolution scenarios discussed in § 4.2. 
The validity of assumption (4) is largely unknown, since there 
is virtually no observational constraint on how much the 
coronal loop area might change between the footpoint and 
summit In the solar case, however, the observed loop area 
variations are small (KJimehuk et al. 1992). We have ignored 
possible area variations and assumed a uniform cross-sectional 
area. We note tins as an uncertainty in our model calculations. 
Assumption (5) may not be satisfi^ in the earliest stages^’of a 
flare, especial^ if the heating rate is very high and the initial 
corona] density is low (Fisher 1986). However, after one or two 
loop acoustic transit times, and after some of the initial evapo- 
ration has taken place, the increased coronal density will 
ensure that the heat flux approaches that given by the Spitzer 
formulation. The long rise times of our flares when compared 
to the loop acoustic transit time justify using the Spitzer for- 
mulation in our models. Assumption (6) is based on the high 
efficiency with which thermal conduction acts to smear out 
coronal temperature variations brought about by spatial 
variations in the flare heating rate. Further, Craig, McCly- 
mont, & Underwood (1978) have shown that the global 
properties of coronal loops (the pressure, apex temperature, 
etc.) are sensitive mainly to the total energy flux deposited in 
the corona and are insensitive to the placement of the energy 
deposition in the corona. Assumption (7) is reasonable for 
coronal loops at temperatures observed during flares, T > 10^ 
K(RTV). 

With these assumptions, the flare loop energy equation for 
A the time rate of change of the loop pressure P, becomes 


3 a ^ „ 5 „ dr dF, 


(4-1) 


— icq is the conductive heat flux and Kq 

(<^ 10'^ in cgs units) is the Spitzer coefficient, v is the plasma 
velocity in the enthalpy flux term, Q is the volumetric flare 
heating rate, and R is the optically thin cooling rate, assumed 
to be of the form R = an^T*. We take a=-^ and 
a = 1.2 X 10~*^ (in cgs units) as reasonable approximations 
(see RTV) to more detailed calculations of coronal cooling 
rates (e.g., Raymond & Smith 1977; Landini & Monsignori- 
Fossi 1990) for these temperatures. It is important to note that 
we have not conadered the effects of possible differences of the 
elemental abundances (especially Fe) from solar values. If it is 
determined that, e.g., Fe has a significantly different abundance 
on AD Leo than on the Sun, this could change the behavior of 
the plasma cooling function and the values of a and a adopted 
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here. In addition, there is stnmg evidence that the quiescent 
coronal abundances of Fe (and other low hrst-ionization- 
potential elements) are enhanced from their photospheric 
values (see, e.g., Meyer 1991). It is not yet clear whether the 
flaring plasma in compact flaic loops on the Sun shares this 
enhancement. In light of these complications, we have simply 
adopted solar photospheric abundances throughout our 
analysis. 

The distance z is measured along the loop from the chromo- 
spheric base, where the temperature T assumes a chromo- 
spheric value Tc^ ; at the loop apex, a distance L (i.e., the loop 
length ”) away, the temperature attains its maximum value 7^, 
with ^ T'lic conductive flux vanishes at both z = 0 

(the chromosphere) and z ^ L (the loop apex). We adopt an 
equation of state appropriate for fully ionized hydrogen, 
P = 2nkT, which relates the number density n of hydrogen 
nuclei to T and P (k is the Boltzmann constant). To find the 
velocity v, we assume that the logarithmic rate of change of the 
density (n/n) along the loop is u^orm and equal to the evapo- 
ration (or condensation) rate N divided by the total column 
number density V, of the loop. This is equivalent to assuming 
that the loop fills or drains homologously. The continuity 
equation can then be solved analytically for the velocity, 
resulting in 

v(z) = N/n{z) X (1 - N{z)/N,) , (4.2) 

where N(z) = n{z')dz\ and N, = N{z = L).^ 

If equation (4.1) is integrated from z = 0 to z = L, the result 
is (see FH90) 

iP = Q-R, (4.3) 

where Q and R denote loop-integrated spatial averages of Q 
and R {Q = Qy since Q is assumed uniform). There are two 
limiting types of behavior, depending on the relative ampli- 
tudes of Q and £ : If 0 |> R, as one might expect early in a flare, 
then radiative losses are unimpoftant in determining the evolu- 
tion of the loop, and the loop nmfergoes ^ strong evaporation,*' 
with the pressure in the loop iocrea^ng as the time integral of 
Q, On the other hand, if £ |> ^ as might happen late in a flare 
after much plasma has been evaporated and the heating rate 
has substantially decreased, then the loop pressure decreases in 
proportion to the average radiative cooling rate, resulting in 
strong condensation.** We modd the rise phase of flares in the 
strong evaporation limit and die decay phase in the strong 
condensation limit. 

Note that if £ = Q, the loop is in energetic equilibrium and 
is therefore a “ static** loop, obeying the loop scaling laws rela- 
ting apex temperature, pressure, and loop length (RTV; Craig 
et al. 1978). The thermal structure and X-ray emissivity of 
static loops are discussed in detail in HF92 and in Hawley & 
Fisher (1994). Static loop modds are used to describe the flare 


^ Similar, but simplified, version 4# CMtr energy equation, applied to the 
decay phase of flares on II Peg and Ainoi were described by van den Oord & 
Mewe (1989, hereafter VM) and Doy^ van den Oord, & Kellelt (1992, here- 
after DVK). In VM, an energy equatkn similar to eq. (4.1) was used, except 
that the enthalpy flux divergence tem was dropped, in which case the energy 
equation results in exactly the same Ikennodynamic structure as that for a 
static loop. On the other hand. VM did account for variations in loop area, 
which our treatment neglects. In OVIC, an order-of-magnitude solution to the 
continuity equation, plus an assumptiOBof exponential decay of thermodyna- 
mic variables with time, were used to define a constant velocity'* cooling 
model for flare decay, which could then be compared with the observed decay 
of temperature and emission measure 


plasma at times between the rise and decay phases of flare EFl 
(Le^ during the flare peak; see § 4.2). 

4,1.1. The Rise Phase 

The loss rate £ is neglected during the rise phase, so that if 
evaporation continues throughout the rise time t,., integration 
of equation (4.3) yields a pressure given by (3/2)£^^ = Jo^ 

Q{t)dt, Here we identify with the elapsed time between the 
onset of the flare rise seen in the EUVE/DS count rate and the 
time of the peak count rate. Because the pressure at elapsed 
time depends only on the time integral of Q, we can replace 
the integral on the right-hand side by Q,, t,., where is the 
average value of Q over the flare rise. Thus, the total energy 
deposited (per unit volume) in the loop (r^Q^y) is equal to the 
volumetric thermal energy at the end of rise phase, (3/2)P^„. 

We can now make use of analytical solutions for the tem- 
perature and density structure in a strongly evaporating loop 
(from Appendix A of FH90) to eliminate Q^y in terms of the 
apex temperature and the loop semi-length L: Q.v = 
where C = 0.2913. Results from Appendix A 
also allow us to eliminate in terms of L, and (the 
column depth in one leg of the loop at the end of the rise 
phase), so that 

Pn^ = 2^'^r,2kT^N^JL, (4.4) 

where q = 0.4656. After we make these substitutions, the solu- 
tion to the time-integrated energy equation r^Q^y ^ (3I2)P^^^ 
reduces to 


_ C>Cq T^a\ 
2^'^n2kL 


(4.5) 


which describes the amount of chromospheric plasma evapo- 
rated into the corona during the rise phase of the flare. 

4.1.2. The Decay Phase 

In the strong condensation limit, appropriate for tl^e decay 
phase of a flare, equation (4.3) reduces to (3/2)P = — £. In § 3 
of FH90, we showed that a good approximation to the spatial 
average of the loss rate is 

R = kN^Jl} , (4.6) 

whwere R.II} is identical to the quantity K defined in equation 
(13)of FH90, Nj is the total coronal column depth in the loop, 
and A = 14/(7 — 2a). Recall that a = —0.5 and thus A = 1.75. 
The value of ^ depends only on the physical constants a, ol, k, 
and defined in § 4.1.1. We relate the condensation rate 
through the loop footpoint N, to P with the Ansatz 

N/N, = QP/P . (4.7) 


Results from a series of numerical gasdynamic simulations of 
flare loop decay (Jakimiec et al. 1992) suggest that Q ^ 0.3 
during flare loop condensation. An interesting discussion of the 
“correct** value of Q during flare loop decay can be found in 
Cargill, Mariska, & Antiochos (1995). Adopting Q = 0.3 and 
combining equations (4.3), (4.6), and (4.7) gives the single differ- 
ential equation, 

dx (2/3)n/eN;i.„ 

dt 

Here X = and and P„„ are the values of the 

column depth and pressure prior to the onset of strong con- 
densation. If strong condensation occurs for a decay time 
after the loop has attained the column depth and pressure 
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and P^, integration of equation (4.8) then yields 

(2/3)flJgiVL.T, 

1/n-A ’ 

where is the value of x after a time has elapsed. 

We now seek to express in terms of an observed quantity, 
in this case the EVVE/DS count rate. In principle, the 
observed EUV emission depends on both the temperature and 
emission measure (hence column depth or density) of the loop 
pia^sma To Understand these dependencies, we have computed 
theoretical EUVE/DS count rates for a sequence of detailed 
loop models in the strong condensation limit (see Appendix B 
of FH90), with each loop model having a different apex tem- 
perature in the range from 1—3 x 10^ K, but with all models 
containing the same coronal column depth. We use the EUV 
plattma emissivities of Landini & Monsignori-Fossi (1990) and 
the EUVE/DS detector response curve from CFAS. The emis- 
sion in the EUVE/DS passband is roughly the same for the 
entire sequence of models, indicating that it is relatively insen- 
sitive to apex temperature variations within this range. The 
emission scales as Nf, or equivalently, x^. Therefore we can 
write where is the EUVE/DS count rate 

at the flare peak, and Cj is the count rate after a time r^ has 
elapsed. 

4 . 1 . 3 . Loop Lengths 

By approximating and in equation (4.9) with 

values of P and achieved at the end of the rise phase, 

equations (4.5) and (4.9) can be combined to yield a single 
equation for the loop length L, which involves only Xj, r^, r„ 
physical constants, and (as estimated from the EUV spec- 
troscopic data): 

1 r (2/3)^nT^Cxo T^^\f-\\/Q - 

^ "" L J 

(4.10) 

Ignoring the .Xj dependence for the moment, and substituting 
A = 7/4 into the expression, equation (4.10) shows that Loc 
^3/7^4/7-pW2 -pjjyj j|,g loop length is relatively insensitive to 
f^, so that modest uncertainties in the apex temperature of the 
loop have a small effect on derived values of L. ^uation (4.10) 
also shows that long rise and decay times require large loop 
lengths. Once L has been determined from equation (4.10), 
estimates for and P^ can be obtained from equations 
(4.5) and (4.4). We emphasize that none of the quantities 
appearing in equation (4.10) require knowledge of the absolute 
flux in the DS detector; relative values of the count rate are 
sufficient to derive estimates of the loop length, coronal 
column depth, and loop pressure. 

4.2. Coronal Evolution During the EUV Flares on AD Leo 

From the EUVE/DS light curve for flare EFl we obtain a 
rise time r, ~ 1 1 ks (see Figs. 1 and 2c and Table 2). TTie 
quiescent count rate observed before the flare is ~0.27 s‘‘, 
while the count rate near the flare peak is 0.98 s ' This implira 
a maximum flare-induced count rate of ~ 0.71 s' ‘. Flare emis- 
sion remains elevated with slight variations for nearly 4 hr 
since the count rate during the peak is nearly constant, we 
assume that the coronal pressure and column depth maintain 
their maximum values during that time. The flare then decays 
rapidly for a time ~ 6 ks. After this, there is a continuing 
slow decay for most of the next day, but our models are only 


apiflicaUe during the initial rapid decay phase. At the end of 
that ph aic, the count rate has fallen to 54% of its p>eak value, 
thus, xj ~ 0 l 54. From the EUV spectrum (Cully et 1995), we 

estimate the peak coronal temperature during flare EFl to be 
7^ ~ 2.5 X 10^ K. Given these parameters, we obtain from 
equation (4.10) a loop length L ~ 3.8 x 10‘® cm (note this is 
frim equation (4.5) we obtain a peak coronal column 
depth of JV,„ ~ 1.3 x 10^* cm'*; and from equation (4.4) we 
obtain a peak cowmal pressure of P ~ 180dynecm'*. 

The quantities L, N and P were estimated using only 

rdatne EUVE/DS count rates. If an absolute count rate is 
available; the loop cross-sectional area and hence the flare 
volume K, the coronal emission measure EM, and the peak 
coronal thermal energy E,h can also be determined by compar- 
ing the observed flux with that found from the loop models. 
We have cmnputed the EUVE/DS flux expected at Earth 
during the flare peak using a static coronal loop with the 
values of L and determined above, coronal emissivities 
taken from landini & Monsignori-Fossi (1990), the DS detec- 
tor re sp o n se curve taken from Figure 2 of CFAS, and estimates 
of the ISM opacity from Rumph, Bowyer, & Vennes (1994) for 
an interstellar hydrogen column depth of 1 x 10^^ cm~^ 
(Paresoe 1984), and a distance to AD Leo of 4.9 pc (Gliese & 
Jahretss 199U The model calculations assume that none of the 
flare pla«na radiation is absorbed by the star. We assume that 
the static loop model applies during the flare peak, siij^c^ the 
count rate remains nearly constant for several hours. The 
average cross-sectional area of the loop is then simply the ratio 
of the observed count rate (counts s”^ and the computed 
count rate (counts s ” ^) for a loop of unit cross-sectional area. 

For flare EFl, the loop cross-sectional area A = 9A x 10^^ 
cm^ which amounts to 0.85% of the stellar surface area if we 
adopt =* 19 X 10‘® cm for AD Leo. (This radius was deter- 
mined first matching the measured spectral type of AD Leo 
to the corresponding effective temperature, then using the bol- 
ometric luminosity and effective temperature to find as in 
Pettersen 1980). The resultant flare loop volume is 
V = 2LA = 7.1 X 10^® cm^, the peak emission measure is 
EM,^ =* {N^L)^V = 8.2 X 10*^ cm'^ and the peak 
thermal energy Eyi - (3/2)P^K = 1.9 x 10^^ ergs. For EFl, 
E,|| is roughly 5 times the observed EUVE/DS energy (E^s in 
Tabic 2). Eih represents the total thermal energy contained in 
the evaporated plasma at the flare peak. In our model, the 
energy radiated during the strong condensation phase orig- 
inates from this thermal reservoir. 

Flare EF2 (Figs. 1 and 2e, Table 2) was both less intense and 
of shorter duration than flare EFl. The “quiescent” level for 
EF2 indodes some post flare emission from EFl, resulting in 
an elevated preflare count rate of 0.46 s From the EUV 
spectrum of flare EF2 (Cully et al. 1995), we estimate a peak 
coronal ^lex temperature « 2 x 10’ K.. The peak flare 
count rate was roughly 0.76 s“S resulting in a flare-induced 
count rate of 0.30 s' *. The rise time for EF2 is poorly deter- 
mined, we can provide only an upper limit, since the count 
rate reaches its peak value only one orbital data gap after being 
at the quiescent level, i.c., ;S 6 ks. The decay timescale is also 
poorly determined for the same reason, namely, the count rate 
returns to its preflare value one orbital gap after being near its 
maximum. Thus wc can set only an upper limit < 6 ks. 
Because there are no reliable measurements for the count rate 
during the decay phase, wc can say only that xj ^ 0 after 6 ks. 
With generous allowance for the uncertainties, we use these 
limits to estimate L, and from our models. We find 
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L - 1.5 X 10'° cm, - LO x cm*•^ and - 280 
dyne cm ~ The estimates arc essentially unconstrained but are 
physically reasonable. Matching a detailed static loop model 
with these parameters to a fiaie-induced count rate of 0.30 s~ ' 
yields a loop cross-sectional area A = 1.9 x 10'® cm^ (or 
0.18% of the stellar sur£scc area), a flare volume 
F = 5.6 X 10^® cm^, a coronal emission measure 
EM = 2.5 X 10^' cm“^, and a peak thermal energy £,h = 2.4 
X 10^^ ergs, or about 4 times tte observed EUVE/DS energy 
in Table 2. 

The derived loop parameters for flares EFl and EF2 are 
summarized in Table 3. 

4.3. Coronal Evolution of the AU Mic Flare 

We can apply the same analysis used to determine the AD 
Leo flaring coronal parameters to the flare peak behavior in 
the giant AU Mic flare observed by EWE and reported in 
Cully et al. (1993). The distance to AU Mic is 9.3 pc (Gliese & 
Jahreiss 1991) and we estimate the stellar radius, using the 
same method as for AD Leo, to be 5.5 x 10'° cm. The 
interstellar column depth of2 x 10'* cm~^ is again taken from 
Parescc (1984), from which the EUV extinction is computed 
using Rumph et al. (1994). Cully et al. (1993) estimated the 
temperature near flare peak to be 7^ ~ 3 x 10^ K. Table 2 
gives our values for the measured rise and decay timescales and 
the corresponding count rates. Note that our analysis applies 
only to the rapid rise and decay seen during the first large peak 
between 197.46 days-197.60 days (sec Fig. 2 of Cully et al. 
1993). Using the published count rates, we find a coronal loop 
length L = 2.6 X 10'° cm, a peak pressure ~ 350 dyne 
cm"^, and a maximum column depth of = 1.5 x 10^' 
cm"^. Using the observed count rat^ we find a coronal loop 
cross-sectional area of 4.5% of the stellar surface area, which 
results in a coronal area coveraw of 1.7 x 10^' cm^, an emis- 
sion measure of 2.9 x 10^^ cm * , and a peak thermal energy of 
4.6 X 10^^ ergs. These values are included in Table 3. Evidently 
the AU Mic flare was much more energetic than the largest AD 
Leo flare (EFl) although the flaring loop (or loop arcade) was 
more compact. We also found that for AU Mic, a static loop 
model with the derived loop length had too large a loop pres- 
sure compared to the value found from equation (4.4). A 
“ strongly evaporating ” loop, on the other hand, gives a consis- 
tent pressure and was used to determine the quoted cross- 
sectional area and volume. A strongly evaporating loop is also 
consistent with the rapid rise to the peak count rate during the 
AU Mic flare, as compared to flare EFl (which rose more 
slowly and stayed at the flare peak value for considerably 
longer, and was best fit by a static loop)L 

4.4. Loop Sizes and Flare Geometry 

The coronal loop lengths and areas we find for flares EFl 
and EF2 (Table 3) are considerably larger than those deter- 
mined for large solar flares, in which a flare loop might have a 
length of a few x 10® cm and a total area of less than 10"^ of 
the solar surface. This qualitative difference between solar and 
stellar flare loops is well known, although our loop lengths are 
somewhat longer than estimates from earlier studies (L 1- 
2 X 10'° cm). Because our technique for estimating loop 
lengths and cross-sectional areas is new, we now compare and 
contrast our methods with those obtained previously. We also 
comment on other interpretations of flare decay that do not 
involve loop condensation. 

Most previous estimates for the sizes of flaring coronal loops 


on dMe stars have been possible only when high-energy 
(usually soft X-ray) data were available, as in Haisch et al. 
(1987), Doyle et al. (1988), and Tagliaferri, Doyle, & Giommi 
(1990). The method used by these authors involves first deter- 
mining the emission measure of the flare from the total mea- 
sured luminosity coupled with an estimate of the emissivity 
(based on an assumed or measured temperature plus an 
adopted cooling law). Next, the mean density in the flare 
plasma is derived by equating the measured flare decay time to 
the radiative cooling rate (again assuming a cooling law), at 
fixed temperature. With the emission measure and density 
established the volume is obtained. To derive a loop length, an 
aspect ratio between the loop diameter and the loop length is 
assumed^ and then the measured flare volume can be related to 
the loop length. Typically, solar values of the aspect ratio (0.1- 
0.2) are invoked, resulting in loop lengths of 1-2 x 10‘°cm. 

Applying this analysis to our data for EFl results in a loop 
length of 4 x 10' ° cm (similar to our model result in Table 3), a 
mean density of 9 x 10'° (about 2.5 times greater than our 
value), and hence a smaller area coverage. The reason for the 
difference in the density is that our condensation model 
accounts sclf-consistently for changes in both the temperature 
and density as the flare evolves. 

In contrast to this earlier work, our estimate of the loop 
length (eq. [4.10]) depends only on measured timescales and 
the coronal temperature and is determined independept}y of the 
volume emission measure. Because of this, the emission 
measure can be used to find the flare volume and the total loop 
cross-sectional area. In our analysis, no assumptions about the 
loop aspect ratio are necessary. Rather, the total cross- 
sectional area is a product of our analysis. If one chooses to 
assume that all loops have a solar aspect ratio of 0.1 -0.2, the 
total cross-sectional area from our analysis can then be used to 
find the number of such loops that comprise a flaring coronal 
arcade. If taken as single loops, the aspect ratios for flares EFl 
and EF2 would be 0.15 and 0.18, respectively. We emphasize 
again that equation (4.10) shows that a flare with a long rise 
and decay time must take place within a long loop or within an 
arcade of long loops. 

Our technique for estimating the loop length is most similar 
to the method described by Reale et al. (1988). These authors 
computed a series of X-ray light curves from several different 
flare loop hydrodynamic simulations, each assuming a differ- 
ing loop length. These light curves were then compared to the 
observ^ light curve from the 1980 August 20 flare observed on 
Proxima Centauri by the Einstein satellite. The energy input 
phase in their simulations was determined by the observed rise 
phase in the flare; the behavior of the decay phase was then 
determined by the assumed loop length. They found that if the 
loop length was too long, the computed decay was slower than 
that observed; if it was too short, the computed decay was too 
fast. This technique for finding L is basically the same as ours, 
except that we have exploited semianalytic models for the flare 
rise and decay phases, rather than employing detailed hydro- 
dynamic simulations. 

Previous work on the active prc-main-sequence dwarf AB 
Dor (Collier-Cameron, & Robinson 1989) suggested that 
centripetal acceleration was required to support the very large 
loops inferred from their data. The loop lengths we find for AD 
Leo arc not as large; coupled with the higher gravity and 
slower rotation speeds of AD Leo and AU Mic, we find that 
centripetal forces are not important for the flare loops in our 
model 
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4.5. Alternatives to Simple Evaporation and 
Condensation in Loops 

If the flare rise and decay phases are a result of different 
physical scenarios than the ones we have assumed in this 
paper, then the interpretation of the EUVEfDS light curve 
evolution will differ accordingly. Two specific alternatives to 
simple loop evaporation and condensation deserve comment, 
namely, reconnection in two-ribbon loop arcades (Poletto, 
Pallavicini, & Kopp 1988, hereafter PPK), and modeling the 
flare decay as the ejection of a coronal loop or plasmoid 
(CFAS). In the former case, PPK fitted the decay phases of two 
stellar flares (the 1985 August 6 flare on EQ Peg observed by 
EXOSAT, and the 1980 August 20 flare on Prox Cen observed 
by the Einstein satellite) with the coronal arcade reconnection 
model of Kopp & Poletto (1984, hereafter KP). In that model, 
magnetic field lines that have been blown open early in the 
flare reconnect into closed field lines. As time proems, a 
sequence of closed coronal loop arcades is formed, in which the 
most recently formed arcade encloses those formed earlier. 
PPK a«inne that the reconnection energy release rate predict- 
ed by KP can be matched in time with the observed X-ray light 
curves. A form for the temporal variation of the reconnection 
rate in the model calculations must be assumed, but this choice 
can be adjusted in such a way as to obtain good agreement 
between the theory and the observations. 

The direct comparison between the theoretical model and 
the observations in PPK glosses over several important issuw. 
Assuming that the energy release rate and the X-ray luminosity 
track in phase with each other is tantamount to assuming that 
the corona can respond instantaneously via evaporation and 
condensation to the heat deposited in the reconnected loops. 
We have shown in § 4.1 that this is not the case. If coronal 
evolution is governed by the KP reconnection model, then the 
response of each one of the coronal loops (or equivalently, each 
arcade) is governed by equations similar to those we discussed 
in § 4.1. Evolution of the entire active r^ion with time would 
be determined by a convolution between the energy release 
rate due to reconnection and the X-ray response function for 
each set of loop arcades. Using a solar example as a guide to 
what might occur during a large stellar flare, consider the 1980 
May 21 two-ribbon flare observed by SMM (Antonucci, 
Gabriel, & Dennis 1984). Even though the decay phase does 
appear to involve the reconnection of open field hnes into 
nested sequences of closed coronal Iooim (see Poletto & Kopp 
1985, 1988), we have found that a single-loop model ade- 
quately describes the evolution of the brightest flare plasma, at 
least through the early part of the decay phase (Fisher & 
Hawley 1990). We believe this is because near the peak of the 
flare, the brightest loop(s) dominate the emission from the 
entire active region, and these loops have only one particular 
length scale. 

An alternative flare decay scenario has been suggested by 
CFAS in their study of the EUVEIDS light curve of the 1992 
July 15 flare observed on the dMe star AU Mic. This flare 
exhibited a brief, intense period of emission (which we dis- 
cussed in § 4.3), lasting roughly 2 hr, followed by a much 
longer, slow decay phase lasting roughly 1-2 days. The two 
different timescales suggest that two different phenomena 
occurred during that flare. For the slow decay phase, CFAS 
propose a model that is based on the ejection of bright coronal 
loops or plasmoids from the star, accompanied by their expan- 
sion into the interstellar medium. While this mechanism could 
explain very long flare decay times, it cannot cope with the 


shorter decay timescales of flares EFl and EF2 on AD Leo (or 
the rapid decay phase of the AU Mic flare), because the 
required expaonon velocities for the ejected loops would be 
unpbysically large, and/or the initial densities would be suffi- 
ciently large that radiative cooling would dominate expansion 
cooling, vioiati^ one of the precepts of their model. Therefore 
wc believe that the initial flare decay seen during the 1992 July 
15 AU Mic flare, as well as the decay phases of flares EFl and 
EF2, ate by cooling and condensation within closed 

coronal loops. 

5. DISCUSSION 

Analysis td the optical data alone allowed us to place the 
observed AD Leo flares in the context of previous optical flare 
data, with the result that they are unusual in terms of longevity 
and relative^ small enhancement but large total energy, 
thou^ they are by no means extraordinary. The optical con- 
tinuum and line emission evolve together in typical fashion, 
and a Uackbody fit to the continuum gives temperatures and 
area coven^es -simila r to other flares. The EUVE/DS data, 
together with our flare models, give an improved determi- 
nation of the coronal loop parameters and flare geometry. By 
combining iilbfmation from the optical and El/FE/DS data 
sets, we now construct a picture of flare EFl as a whole and 
compare it with what is known about solar flares. 

5.1. A Picture Flare EFl 

We adopt a coronal loop geometry with footpoints in the 
photos|diere and height ~ 1 H, above the stellar surface. We 
refer to the “loop" as a single entity but do not exclude an 
arcade of loops of similar length whose total area combines to 
give the inferred area coverage. 

The flare is fint seen as a rapid increase in the optical contin- 
uum, observed in the 17-band photometry. Hudson et al. (1992) 
showed that whHe-light and hard X-rays in solar flares were 
coincidait both in space and time. Since the hard X-rays are 
emitted 1^ a cc ele rated nonthermal electrons, wc identify the 
optical continssm radiation with the impact sites where the 
electrons are depositing their energy via collisions in the dense 
lower atmosphere. Our analysis of the optical emission indi- 
cates that the fractional area coverage (~0.01% of the stellar 
area) is much smaller than the estimated coronal cross- 
sectional area (~ 1%). On the Sun, the chromospheric flare 
kernds ima ged in Ha or hard X-rays are also found to be much 
smaller in extent than the coronal flare loop areas (Canfield et 
al. 1991). One possible explanation for these differences in area 
coverage is the convergence of magnetic field lines in the 
chromosphere and photosphere, which results from the tran- 
sition between a space-filling potential or force-free field in the 
corona, to thm, discrete magnetic flux tubes confined by gas 
pressure at the photosphere. 

Since EFl occurred in a rather long loop, the initial flare 
heating seen in the U band did not produce an immediate 
change in the coronal emission measure, although evaporation 
was probably beginning at this time. 

The second stage of the flare is characterized by a plateau in 
the (/-band emission which we identify with a period of contin- 
ued electron picci[ntation possibly at fresh impact sites within 
the loop. The loop is undergoing strong evaporation and the 
coronal emisskm begins to increase. 

At the end of this period, the (/-band emission again rises 
sharply to a second peak, while the coronal emission reaches 
its peak. The combination of optical and £(/K£/DS data over 
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the entire flare to this time (i.c„ the time of maximum coronal 
emission measure), now allows us to examine the evidence for a 
stellar Neupert effect. As discussed in § 1, the canonical solar 
model leads to the prediction that hard X-rays and white-light 
emission will precede the soft X-ray (EVVEfDS) emission and, 
in fact, the EUVEfDS emission should scale with the time 
integral of the white-light energy. The time delay we observe 
between the optical flare and the coronal emission could there- 
fore be interpreted as a stellar Neupert effect. In Figure 3 we 
show the time-integrated U-band energy superposed on the 
observed £l/F£/DS light curve (with the count rate expressed 
in energy units). The temporal behavior matches well and, if 
anything, shows that the nonthermal electron power (which 
must be at least as great as the power observed in U) is more 
than adequate to produce the observed EUV (i.e., coronal) 
emission during the flare evolution. In § 5.3 we discuss some 
previously observed flares and find that they also appear to 
show a Neupert effect. Evidently, the same underlying heating 
processes, direct excitation by accelerated electrons and evapo- 
ration of chromospheric material into the corona, can account 
for the temporal relationship between coronal and optical con- 
tinuum emission on both dMe stars and the Sun. However, 
stellar flares can be far more energetic and exhibit heating 
timescales that are significantly longer than those observed in 
the solar case. We discuss some possibilities for why these 
differences might arise in § 5.2 below. 

Following the rise phase, the coronal emission reaches a 
plateau near its maximum value for a period of several hours. 
The radiative cooling time in the corona during this period is 
about 10^ s; since the emission shows no decay over that time- 
scale, we surmise that there must be continued heating. The 
heating is not accompanied by evaporation, since the emission 
measure does not increase. One possibility may be that the 
flare heating mechanism which formerly reached the lower 
atmosphere is now depositing its energy directly into the 
corona; we address this further in § 5.2. Unfortunately, we 
have no optical observations to confirm the lack of lower 
atmospheric heating during this phase. 

We note that the poor time resolution of the EUV data 
allows an alternative explanation for the emission during this 
time period. It is possible that the “flare peak” between ~36- 
40 hr is actually a composite of emission from two flares, one of 



Fig. 3. — Solid curve shows a time integral of the (/-band luminosity (with 
(he prepare luminosity subtracted) during the rise phase of flare EFI. Filled 
circles show the EVVE/DS count rate during this penod, converted to lumi- 
nosity units. Wc argue that the behavior shown here provides evidence for a 
stellar Neupert effect. 


which peaks at ^36 hr, starts to decay in the next orbit, and 
then is joined by a second flare which peaks at ^40 hr. This 
would remove the need for a continuous heating process 
lasting some 4 hr at the flare peak. Our data do not allow us to 
distinguish between this two-flare scenario and a single-flare 
plateau phase. However, if we assume a plateau, our analysis of 
the loop during this phase shows that a static model is consis- 
tent with the loop length and pressure derived from equations 
(4.10) and (4.4), which does support a slow deposition of heat 
directly into the corona. 

The coronal emission next enters a phase of relatively rapid 
decay, which we identify as the strong condensation phase. 
Flare heating has ceased and the coronal plasma cools and 
condenses down the loop into the chromosphere. It is the com- 
bination of the rise and decay phase timescales, together with 
our theory of evaporating and condensing loops, which allows 
us to determine the loop length unambiguously. Thus, the long 
loop length we derive follows directly from the observed long 
rise and decay times. As mentioned earlier, this does not pre- 
clude a scenario in which several episodes of heating take place 
within an arcade of loops of this characteristic length. 

It is possible that an ensemble of short loops (which would 
respond more rapidly to flare heating), together with a con- 
trived heating profile, could mimic the natural response of a 
long loop. We regard this as an unnecessary complication 
unless further observational constraints on the loop lerjgth are 
discovered. ' 

Finally, the coronal emission switches from a (relatively) 
rapid decay to an extended slow decay, lasting most of the next 
day. Our models do not address this phase. It may represent a 
relaxation process as the loop returns to quiescence, or, as 
suggested in CFAS for AU Mic, it may be the emission from an 
expanding and cooling coronal mass ejection. One intriguing 
feature of the slow decay phase on both AU Mic and AD Leo 
is the presence of a second small flare (here, EF2) which 
occurred approximately 1 day after the initial large flare (here, 
EFI). This may be merely coincidence, though Brown et al. 
(1994) also report a similar secondary flare on HR 1099 
approximately I day after a large flare. Further observations 
are required to test whether the second flare is a typical (and 
perhaps important) feature of stellar flares. 

5.2. Impulsive Phase Quenching 

As discussed in § 4, the thermal response times and the 
amounts of energy which can be stored in hot coronal plasma 
should depend in a predictable way upon the length of coronal 
loops. Wc showed, in particular, that the loop length L cc 
A more subtle issue is whether the coronal loop length 
and/or cross-sectional area can have a direct effect on the dura- 
tion of the impulsive phase of flares. Such a relationship is 
suggested by the longer duration of impulsive optical contin- 
uum emission in many stellar flares as compared to solar flares. 
One possibility is that impulsive phase timescales are set by an 
intrinsic flare reconnection rate and the amount of magnetic 
flux available to be reconnected. If this were the case, we might 
expect the duration of the impulsive phase to be related to the 
coronal loop cross-sectional area or diameter, which is often 
found to be larger in the stellar case than in the solar case. 
Another possibility, described below, is that the impulsive 
phase of flares ends when it is “ quenched ” by a high coronal 
density. The lime it takes for the quenching to occur is a 
strongly increasing function of loop length. 

The impulsive phase acceleration mechanism, whatever it is, 
must obey at least two constraints: (1) Energy is extracted (by 
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some means) from nonpotential magnetic fields and converted 
into the kinetic energy of the accelerated electrons, and (2) the 
acceleration mechanism must be effective enough to overcome 
Coulomb losses from partidei in the ambient plasma. An 
extended period of flare energy release will result in an increas- 
ing level of the overall coronal density, because of chromo- 
spheric evaporation. If the impulsive phase is to persist, the 
acceleration mechani sm must continue to operate, even in the 
face of the increased Coulomb losses. We suggnt that the 
impulsive phase will be quenched when the ambient coronal 
density reaches a critical threshold where Coulomb losses 
overcome acceleration. The tune it t^es for this to occur 
would then determine the duration of the impulsive phase. 
Because the loop length L affects how quickly corontd loops 
can fill with evaporated plasma, we then expe^ a relationsUp 
between coronal loop length and the duration of the impulsive 
phase. For a given flare apex temperature T^, the quantity 
n = I^,d/L scales as L~^ (eqs. [A8] and [A9] of FH90) in the 
strong evaporation regime. Writing the duration of the impul- 
sive phase gives T,^xL^ meaning that the 

impulsive phase for stellar flares occurring in long loops should 
indeed be much longer than that for solar flares. 

Thus we have identified at least two possibilities for why the 
duration of the impulsive phase should depend on loop length 
or cross-sectional area, and hence why the impulsive phase for 
man y stellar flares might last longer than that for solar flares. 
However, determining which (if either) scenario is operable can 
probably only be addressed by e xa m i nin g existing data sets of 
solar flares, where both good spatial and temporal resolution 
are available. 

The wide range of observed stellar rise times, and the 
varying degree of correlation between optical and high-energy 
li gh t curves that have been cited in previous stellar flares (see 
§ 1), thus finds a natural explanation in terms of a wide varia- 
tion of loop sizes in different stellar flares. 

5.3. Comparison with Previous Flare Observations 

Several previous flare observations show evidence for a 
stellar Neupert effect such as the one reported here. Exami- 
nation of Figures 2 and 3 from Kahler et al. (1982) shows that 
the [/-band emission peaks during a time when the X-ray 
count rate is rising rapidly for two separate flares on YZ CMi. 
A similar result is visible in Figure 11 of dc Jager et al. (1989), 
during a flare on U V Ceti. 

The existence of a physical connection between the coronal 
and chromospheric emission, as manifested by the Neupert 
effect, provides a straightforward explanation for the empirical 
correlations between energy emitted in X-rays, chromospheric 
emission lines, and optical continuum emission discussed in 
Butler et al. (1988), Haisch (1989), and HP91. The energy corre- 
lations extend over several orders of magnitude in total flare 
energy and apply to both the Sun and to flare stars. HP91 
interpreted this result as evidence for a constant heating rate in 
flares of very different total energy, so that the energy difference 
depended on the area coverage and duration of the heating. 
Our new data and models suggest that the primary source for 
energy deposition observed in solar flares, i.e., a flux of non- 
thermal electrons, can also produce the emission seen in stellar 
flares, provided that the flare occurs in a loop of appropriate 
length as determined by the rise and decay timescales. It is then 
the relative efficiency of the nonthermal electrons in heating 
the chromospheric material directly to produce the optical 
emission, and in evaporating it to coronal temperatures to 
produce the soft X-ray emission, that must remain nearly con- 


stant in order to obtain the observed empirical correlation 

b e twee n the ene r gies. 

Evidence for such a constant relative efficiency is indeed 
suggested in the case of solar flares by the work of Kane & 
DoimeDy (1971), Donnelly & Kane (1978), and Kane, Frost, & 
Donnelly (1979X who showed that the ratio of long-wavelength 
EUV radiation (emitted primarily from the transition region 
and diromosphere) to hard X-rays was nearly constant over a 
very wide range of solar energies. A more recent study of white- 
light emission and hard X-rays (Neidig & Kane 1993) showed 
similar resuhaL 

6. CONCLUSIONS 

Out simultaneous optical and EUVE/DS data obtained 
durii^ a flare on AD Leo provide strong evidence for a stellar 
Neupert effect in which the EUVE/DS light curve is roughly 
{voportkmal to a time integral of the l/-band light curve. 
Wba we eiamine earlier simultaneous optical and soft X-ray 
stellar flare data, we find confirmation of this phenomenon. 
The wdl-estaMished spatial and temporal correlation of white- 
light emtssioo and hard X-ray emission in solar flares, together 
with the identification of the hard X-ray emission as a signa- 
ture of nonthermal bremsstrahlung from a flux of energetic 
electrons., suggests that the optical continuum produced in 
stellar flares is a real observational signature of impulsive 
phase nonthermal electron heating. Therefore, the mod^ that 
has be en developed to explain solar flare observafions is 
equafly applicable to stellar flares. Differences in rise and decay 
timescales, and flare energetics, can be ascribed to differences 
in the sizes of the coronal loop structures. 
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